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Abstract 

We give texture analyses of cascade hierarchical mass matrices in su- 
persymmetric 50(10) grand unified theory. We embed cascade mass tex- 
tures of the standard model fermion with right-handed neutrinos into the 
theory, which gives relations among the mass matrices of the fermions. 
The related phenomenologies, such as the lepton fiavor violating pro- 
cesses and leptogenesis, are also investigated in addition to the PMNS 
mixing angles. 
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1 Introduction 



The neutrino oscillation experiments have suggested that there are two large mixing angles 
among three generations in the lepton sector while all mixing angles in the quark sector 
are small. It is known that the current experimental data of leptonic mixing angles [l] is 
well approximated by the tri-bimaximal mixing [2], which is given by 

/ 1/v^ \ 

^TB = -1/v^ 1/v^ -1/V2 . (1.1) 

V -1/V6 1/v^ 1/V2 / 

Such suggestive form of the generation mixing gives us a strong motivation to study a 
flavor structure of the lepton sector. Actually, there are a number of proposals based on a 
flavor symmetry to unravel it and related phenomenologies have been elaborated [3]. 

It has been pointed out that the neutrino Dirac mass matrix of the cascade form can 
lead to the tri-bimaximal mixing at the leading order in the framework of type I seesaw 
mechanism |1] . The mass matrix of the cascade form is parametrized by 



Mcas ~ (5 A A w, with 151 < lAI < 1, (1.2) 




and V denotes an overall mass scale. We call this kind of hierarchy and the matrix with such 
hierarchy, "cascade hierarchy" and "cascade matrix" , respectively. On the other hand, the 
down quark mass matrix of a different hierarchical form, which is 

/ e' 5' 5' \ 

Mhyb ~ 5' A' A' t/, with \e'\ < \5'\ < |A'| < 1, (1.3) 
\5' X' I J 

can give realistic values of CKM matrix elements. The (1,1) element, e', of this matrix 
is smaller than all other elements but this hierarchical structure is close to the cascade 
form except for e'. We call this type of hierarchy "hybrid cascade (H.C.) hierarchy", and 
the matrix with such a hierarchy "hybrid cascade (H.C.) matrix". The neutrino Dirac 
mass matrix of a cascade form gives nearly tri-bimaximal generation mixing and the down 
quark mass matrix of a H.C. form realizes the CKM structure. The fact gives us a strong 
motivation to comprehensively investigate the quark and lepton. Actually, a proposal to 
embed such cascade textures into a supersymmetric (SUSY) SU{5) grand unified theory 
and texture analyses have been presented [3]. For comparison with a SU{5) case and its 
results, we investigate embedding cascade hierarchies into a SUSY 5*0(10) GUT in this 
paper, which is also one of fascinating grand unified models. 

The paper is organized as follows: In section 2, we give a brief review of cascade 
hierarchies for the fermion masses and mixing angles. In section 3, we embed the cascade 
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hierarchies into the SUSY 5*0(10) GUT. The texture analyses for the quark/lepton sectors 
are also given in the section. In section 4, we show some numerical analyses of our model. 
In section 5, we give a comment on the proton decay. Section 6 is devoted to the summary. 
Appendix A gives a discussion about constraints on the structure of right-handed neutrino 
mass matrix. 

2 Cascade hierarchies for fermion mass matrices 

In this section, we give a brief review of cascade hierarchies for mass matrices of the 
fermions. First we discuss the cascade textures for quark and lepton sectors independently. 
The study of cascades for the lepton sector has been discussed in Then a possible 
expansion of the study to quark sector was presented in [5], which was considered in a 
SUSY SU (5) GUT. The work [1] has pointed out that the neutrino Dirac mass matrix of 
a cascade form can lead to the tri-bimaximal mixing at the leading order in the framework 
of type I seesaw mechanism. Since the tri-bimaximal structure can be almost induced from 
the neutrino sector, the mixing angles from the charged lepton sectors should be small [5]. 
This means that the form of charged lepton mass matrix can be taken as a cascade or 
H.C. because mixing angles for those textures are small enough. How about the quark 
sector? The CKM structure is almost determined by a structure of down quark mass 
matrix because of large mass hierarchies among up-type quarks. And it is known that 
the down quark mass matrix of a H.C. form can give the realistic CKM structure (e.g. 
see in E]). The contributions from up-quark sector to the CKM mixing are automatically 
tiny. This means that the form of up quark mass matrix can be taken as a cascade or 
H.C Finally, we comment on the structure of right-handed Majorana mass matrix. The 
contribution from the right-handed Majorana mass matrix should be also small because a 
nearly tri-bimaximal mixing are almost induced from the neutrino Dirac mass matrix with 
the seesaw mechanism, which means it is possible to take the right-handed Majorana mass 
matrix as a cascade or H.C More detailed explanations about the above points including 
mass eigenvalues given from each mass matrix of cascade and H.C. is given in |5]. Here, we 
summarize the above discussions about possible structures of mass matrices of the fermions 
as. 



Mr 



u 



cascade. 



cascade or H.C. or small mixing matrix, 
cascade or H.C. or small mixing matrix. 



cascade or H.C. or small mixing matrix, 
H.C, 



(2.1) 
(2.2) 
(2.3) 
(2.4) 
(2.5) 
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where Mu, M^, M,^/), Mg, and Mr are mass matrices of up- and down-type quarks, neutrino 
Dirac, charged lepton, and right-handed neutrino, respectively. 

Next, we comment on the cascade textures in a SU{5) case. The SU{5) GUT predicts 
a relation between mass matrices for the down-type quark and charged lepton. 

Me ^ Mj, (2.6) 

due to an unification of matter contents. As discussed above, since only mass matrix of a 
H.C form are allowed for Md in the study of cascade texture, the mass matrix for charged 
lepton should also have the H.C. from. On the other hand, some hierarchical structure of 
the mass matrices for the up-type quark and right-handed neutrino are allowed as long as 
induced mixing angles from these matrices can be treated as collections for the CKM and 
PMNS structures, respectively. Therefore, we can parametrize the mass matrices of the 
cascade or H.C. form for the fermions as 

^ I 5. A„ A„ with I !^"! = lf"l^^l^"I^^V.'S'n'''^'' (2-7) 
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Xd 


Xd 


Sd 


Xd 
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5u 




1 



|e„| < \5u\ < |A„| < 1 : H.C, 
Md ~ \ 5d \d Arf I ^dVd, with \ed\ < \6d\ < |Ad| < 1 : H.C, (2.^ 



M^D — \ K \ iuVu, with \5y \ <^ \X„\ <t: 1 : cascade, (2.9) 

Su Xy 

^d Sd 

Me ~ \ 6d -3\d \d I ^dVd, with le^l < \6d\ < \\d\ < 1 : H.C, (2.10) 
Sd Aq 

without 0{1) coefficients for all elements. Here Vu and Vd are vacuum expectation values 
(VEVs) of up- and down-type Higgs fields in a supersymmetric scenario, and the overall 
factor C,d and could be small. We also notice that the Georgi-Jarlskog (GJ) factor [6] 
is introduced to mass ratio between the down-type quarks and charged leptons for each 
generation, 

"^^^ - ^,3,iV (2.11) 



mt, TTis md / V 3 

3 Cascade hierarchies in 5*0 (10) GUT 
3.1 A SUSY 50(10) Model 

We consider about embedding the (hybrid) cascade hierarchical mass matrices into 5*0(10) 
GUT in this paper. A simple 5*0(10) GUT predicts relations between mass matrices the 
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up-type quark and neutrino Dirac, 

M„ ^ M,D, (3.1) 

in addition to the relation (12. 6p . As discussed above, since only the mass matrix of the 
H.C form is allowed for M^, the mass matrix for the charged lepton should also have the 
H.C form like in SU{b) case. For the up-type quark sector, a simple 5*0(10) case predicts 
a GUT relation of the mass matrices ~ M^o-, and the up-type quark mass matrix 
should be restricted to a cascade form because the cascade form of neutrino Dirac mass 
matrix is needed for generating the tri-bimaximal neutrino mixing at the leading order. For 
the structures of right-handed neutrino mass matrix, some arbitrary matrices are allowed 
as long as induced mixing angles can be treated as collections for the PMNS matrix. 

To demonstrate the idea, we consider a simple SUSY S'O(IO) model, which the Stan- 
dard Model (SM) fermions with the right-handed neutrino are included into the spinor 
16- dimensional representation, ip. To give suitable fermion masses, we introduce the Higgs 
fields, i.e. two Higgs 10-plets, H12 and two Higgs 126-plets, Ai^2- There are several ways 
to break 5*0(10) down to the SM. Here, we consider a minimal framework where the break- 
ing of 50(10) is achieved by the Higgs 210-plet [3, [8], $, which breaks 50(10) down to 
Pati-Salam group: SU{A)c x SU{2)l x SU{2)r. We choose that the Pati-Salam group 
is broken further down to the SM via the VEV of the SM singlet component in A2 and 
the VEV also gives Majorana masses for the right-handed neutrinos. Since this singlet 
VEV gives the non-vanishing contribution to D-term in the superpotential resulting in the 
unwanted source of SUSY breaking at high energy (close to the GUT scale), we intro- 
duce a 126-plet, A, whose SM singlet component obtains the VEV to cancel the D-term 
contribution (for instance, see [9l HO]). Here we include two 10-plets, Hi 2 because the 
mass matrices of the up-type quark and down-type quark have to be different in order to 
predict the correct CKM mixing angles, as well motivated from the previous discussion. 
Moreover, we also need one 126-plet, Ai, in order to achieve the GJ relations (12. lip , that 
is, to give the factor of —3 in the (2,2), (2,3), and (3,2) elements of the charged lepton 
mass matrix with respect to that of the down quark mass matrix *] Another 126-plet, A2, 
is introduced to generate the different texture for the right-handed neutrino masses and 
also break the Pati-Salam group to the SM. In our setup, there are six pairs of Higgs dou- 
blets, 0„ = {Hi-a, H2,u, Ai,„, A2,u, A„, and (pd = H24, Ai,^, A24, A^, with 
the mass term (j)uMH(p'd- Note the label u, d refer to the SU{2)l doublet component with 
hypercharge ±1/2 within the GUT multiplet. The mass matrix Mh can be diagonalized by 
UJ^MhU^^, which U^^, U^^ are unitaty matrices acting on 0„ and 0^ respectively. In the di- 
agonal basis, the Higgs fields are given by (0'Ja = (t/^J/3a(0n)/3 and {(p'^a = {U^Jpa{(pd)p- 



* We note that the factor of —3 can be obtained by the couphng of the Higgs 120 or 126, see for 
instance [llj . 
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For the sake of the study, we will not specify how 50(10) is broken in detail, but by 
some doublet-triplet splitting mechanism (for instance see [IS [131 [H]) we will assume that 
Hu = (0m) 1 and Hd = {4>'fi)i have mass at the electroweak scale while the others are so 
heavy and decoupled from the low energy theory. The Higgs fields, a are the two Higgs 
doublets of the Minimal Supersymmetric Standard Model (MSSM). 
The superpotential of the model is given by 

Wy = Y-^'^ipHiiP + f2^Vi/2^ + Y^^^ipAiifj + Y^^^i;A2ilj , (3.2) 

which can be written in terms of the SM components as follows [T5] : 

+L{Yl'Hi,u + Y,''H2,u - SY.'^^Ai^u - 3Yi^''A2,u)N 

+Q{Y,''^H,,, + Yi'H2,d + r/'^Ai,, + F2'''A2,rf)D^ 

+L(Yii°ifi,rf + Y2''H2,d - 3F/^^Ai,, - SYi^'' A2,d)E' , (3.3) 

where the doublet component in the GUT multiplet can be written in term of the MSSM 
Higgs doublets as (0„)o = (U^JaiHu and (0^)^ = (U^JaiHd- 

For the neutrino sector, we assume that the SU{2)l triplet component, A2,t, and the 
SM singlet component, A2,S; in A2, give tiny Majorana masses for the left-handed neutrinos 
and the heavy Majorana masses for the right-handed neutrinos respectively. This results 
in the seesaw formula as follow: 

= Mll - Ml^M^^M.D , (3.4) 

where Mll = Yi~^^{A2,T) = ^^'^^vl, Mr = Yi~^^{A2,s) = Y^^vr and M^d is the Dirac 
mass term whose structure will be discussed below. Since the triplet VEV (A2,t) = vl 
depends on parameters in Higgs superpotential (for instance, see [9]), we assume that the 
VEV is tiny such that the second term in Eq (13. 4p dominates, resulting in the type I seesaw 
dominance. Note that the singlet VEV, (A2,5) = vr, is of order 10^^ GeV. 

After the electroweak symmetry is broken via the doublet VEVs, {Hu,d) = Vu,d, the 
fermion masses are given by 

M„ ~ (f/^Jiin%. = (3.5) 
M,, ^ (f/^JiiF/V = On (3.6) 
Md ^ ((f/^j2iV'2^° + (t/0j3i>'r)^. = (>^2'° + >^r)^'. (3.7) 

Me ^ ((t/^j2ll^2'°-3(f/<^j3in'><^=(>^2'°-3r/>,, (3.8) 

where we assume that the main contribution for the up-type quark (Dirac neutrino) masses 
comes from the coupling to Hi while for the down-type quark (charged lepton) masses 



5 



they arise from the H2 and Ai couphngs. These can be achieved through the following 
assumptions: > (f/</,„)ii and (t/</,rf)2i, (t^^Jsi > (f^^Jii- The Yukawa couplings are 

defined as Y^"" = (f/^Jn^/o, y^^o = (f/^Jai^a'" and F/^e ^ (t/^j3^1>i26_ j^pogg ^^e 
hierarchical forms to the Yukawa couplings, 

Su K K , with \6u\ < |A„| < 1, (3.9) 
Su Xu 1 / 

/ 5d 5d \ 

r^io ~ { 5d 5d 5d \, (3.10) 
^ I A, A, I , with le^l < < \\d\ < 1. (3.11) 



The structure of ¥2"^^ will be discussed in the following sections in term of the right-handed 



neutrino mass matrix, Mjj = Y2^^vr. These lead to the fermion mass matrices at the GUT 



scale as 



M„ ^ \ Su A„ A, Vu, (3.12) 

/ 5u 6u 5„ \ 

Mi,£) 5u Au A„ Vu, with <^ \Xu\ <^ 1 : cascade, (3.13) 
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A„ 1 


ed 


Sd Sd 


Sd 


Ad Ad 


Sd 


Ad 1 


ed 


Sd 


Sd 


-3Ad 


Sd 


-3Ad 



Ma - \ 6d Ad Ad ^dVd, (3.14) 




M, \ 5d -3Ad -3Ad Cd^d, with |ed| < |5d| < |Ad| < 1 : H.C., (3.15) 
V -3Ad 1 / 

where 0{1) coefficients for all elements have been dropped. 



3.2 Cabibbo fitting of cascade mass matrices 

The cascade hierarchical parameters are determined by observed values. It is naturally 
expected that such hierarchies are originated from a symmetry and/or some dynamics in a 
high energy regime rather than solely determined by the magnitudes of Yukawa couplings. 
Although the origin of the hierarchies is not specified in the analysis, one can estimate and 
study the relative magnitudes of the hierarchies introducing a small parameter in the mass 
matrices. In the following, we choose the Cabibbo angle, sin 6'c ^ A = 0.227, as a fitting 
parameter, and study significant implications of the cascade 50(10) scenario. Then we 
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have 



A„ ~ 0.87 X A^ Su ^ 0.85 x A* 



(3.16) 



for up-quark mass matrix of the cascade form and 

\d ~ 0.35 X A^ 6d ~ 0.35 x A^ 



(3.17) 



for down-quark one of the H.C. form at GUT scale, where we utihzed values of quark 
masses listed in [T6] . 

Notice that the C,d is a parameter, which determines a ratio between (3,3) element of 
Yukawa matrices for up- and down-type quarks, and thus, it is correlated with the tan /3 
as. 



tan /3 



Vd 



mt/rrib ~ (9(50) for C^d ~ [large] 
Xmt/mb ~ 0(10) for C^d ~ [moderate] 

small] 



(3.18) 



AV/m,~0(l) for ~ A^ 
As the results we can write cascading textures at GUT scale as 



Md 



X' X' X^ 

X' X^ X^ 

X' X^ 1 

A3 A2 A2 

A3 A^ 1 

A^ A3 A3 

A^ A3 A 

A^ A^ A^ 

A^ A^ A^ 



(3.19) 



Vd [large tan/3] 

Vd [moderate tan /3] 

Vd [small tan (3] 



(3.20) 



where /c^ > 1 is needed to obtain suitable mass eigenvalues after diagonalizing these matri- 
ces. It should be remembered that ~ Md but the additional GJ factor —3 is multiplied 
to the (2,2), (2,3), and (3,2) elements of Me as discussed in the previous section. 

3.3 Neutrino sector 

Next, we consider the structure of neutrino mass matrices. In the cascade model [U [S], 
cascade parameters are constrained as 



5. 



A. 



^""'i ^ 3.19 X 10-2 < A^ 



(3.21) 
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in order to preserve the tri-bimaximal mixing at the leading order with 

Am^i = (7.695 ± 0.645) X 10"^ eV^ (3.22) 
\Aml^\ = 2AOtoil X 10"^ eVl (3.23) 

at the 3a level [1]. Due to the SO{10) GUT relation M^£, ~ M^, the neutrino mass matrix 
can be parametrized as 





A8 


A« 


A« 


A^ 


-A^ 




-A^ 


1 



(3.24) 



where we note that an opposite sign between (2,2) and (2,3) elements is experimentally 
required to obtain the tri-bimaximal mixing as commented in [3]|l] 

3.3.1 Diagonal Mr case 

Let us discuss the case of a diagonal Majorana mass matrix of the right-handed neutrinos, 
Mn = Diag[A^^, A^^, 1]M, where xi > X2 > 0. The cascade model requires the normal 
mass hierarchy of light neutrino mass spectrum in order to realize a nearly tri-bimaximal 
mixing The mass eigenvalues can be estimated as 

mi - ^=^1' (3-25) 
7712 ^ (sX'^-'' + ^Y^=m2 + 2mi, (3.26) 

ms - (^2A«-^^ + ^ = m3 + 3m3, (3.27) 

with a leading order corrections of (9(mi). In order to understand the hierarchical structure 
of the mass matrix and the constraints on the cascade parameters, we write down the 
effective neutrino mass matrix as 




We find that if the terms in the first and second lines are leading contributions, the tri- 
bimaximal mixing can be realized at the leading order. In order that the first term in the 



^ Since the Dirac mass matrix and the up-type quark mass matrix are constrained to have the same 
structure, this opposite sign is also imposed on the up- type quark mass matrix (|3.19p . 



8 



second line does not spoil the structures given in the first line, mi <^ m2, ms is required. 
That is the reason why the neutrino mass spectrum in the cascade model should be the 
normal hierarchy. In the case, we can well approximated as 



7712 ^ Y Am^i and m-s ~ y\Amli\. (3.29) 

Now we can obtain the following four constraints on the cascade parameters: (i) The 
neutrino masses should satisfy rrii ^ This means that Xi > 17 for the parameters by 
utilizing fl3.25p and f l3.26p . This constraint leads to small mass of the lightest right-handed 
neutrino as shown later, (ii) In order to be consistent with a experimental data for the 
neutrino mass squared difference as 



m 



2 

21 



m 



2 I 
31 1 



0.18, (3.30) 



one should have a relation among the cascade parameters as Xi — X2 = 7 or 8, where we 
use the fact that A ~ r. (iii) We have a relation among the cascade parameters, light and 
heavy neutrino mass scales, that is, 

M^4=^, (3.31) 



where 7713 ^ y^jAmfj is taken, (iv) The hierarchy m2 ^ msA^ is required in order that 
the second term in the last line of (I3.28P does not spoil the democratic structure in the first 
line. This gives a constraint Xi — X2 > 5. The above four constraints restrict the neutrino 
Dirac mass matrix of the cascade form and the right-handed one of the diagonal form to 
textures presented in Tab. [HandO We find that the minimal model for the neutrino mass 
matrices is described by (xi,X2) = (17, 10) given in Tab. [1] In this case, mass spectrum of 
the right-handed neutrinos is estimated as 



(Ml, M2, M3) ~ (10^ 10'", 10'^) GeV. (3.32) 

Here we comment on the predicted mixing angles from cascade model. The mixing 
angles of the cascade model deviate from the exact tri-bimaximal mixing angles even if 
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A^ -A^ 1 J 
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( A« A« A« \ 
A8 A^ -A^ 
A8 -A^ 1 ) 






/ Ai« ^ 
All 

^ 1 y 













Table 1: The textures of the neutrino Dirac mass matrix of the cascade form and the right- 
handed neutrino Majorana one of the diagonal from constrained by the experimentally 
observed values of the neutrino masses with the condition x^ — X2 = 7. 

the right-handed neutrino mass matrix is diagonal. The mixing angles can be estimated as 

o 2 



sin^^ 



12 



1 



^12 



sin^ 



23 



1 2 mi 
1 1 



V2 



(1) 

13 



^23 



1 mi(3m3 — 7712) \^ 



1712 



13 



V2 m3{ms 
2 



7712 



3\/2m3 - m2 



-^^13 



in a perturbative methodjil where parameters 91^' indicate deviations from the exact tri 
bimaximal mixing angles. These are elements of the following mixing matrix, 

/ 1 



- 1^2 



771x7712 



y/2m3-m2 m3(m3-m2) 



(1) 



(3.33) 
(3.34) 
(3.35) 
(3.36) 
(3.37) 

(3.38) 



y(i) 



'12 
3(1) 



'^12 
1 



'^13 > 



e. 



(1) 
23 

1 



(3.39) 



In our notation, the experimentally observed PMNS mixing matrix is given by Vrmns — 
Vt:bV^^^Pm, where the Pm is a diagonal phase matrix. 

3.3.2 Non-diagonal Mr case 

We discuss the case of non-diagonal M/j, which is generically allowed in the context of 
the cascade textures. First, we define the diagonalized right-handed neutrino mass matrix. 



^See [5] for a detailed derivation. 
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Xi 


X2 




MrIM 


17 
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1 A8 A8 A8 \ 

A^ -A^ 
\ A^ -A^ 1 / 


/ n n \ 

/A U U \ 

A9 
V 1 / 


18 


10 


/ A« A« A« \ 

A8 A^ -A^ 
\ A8 -A^ 1 / 


/ Ai« \ 

Aio 

V 1 y 











Table 2: The textures of the neutrino Dirac mass matrix of the cascade form and the right- 
handed neutrino Majorana one of the diagonal from constrained by the experimentally 
observed values of the neutrino masses with the condition x\ — X2 = 8. 

Dr, as 

/ A^i \ 

Dr = UIrMrU.r =\ A"2 M with xi > X2 > 0, (3.40) 



where Mr is a non-diagonal mass matrix for the right-handed neutrinos but mixing angles 
among each generation are assumed to be small in order to preserve the tri-bimaximal 
mixing. If the mixing angles among each generation of the right-handed neutrino are small 
enough, U^r can be written by 

/ 1 Qr,i2 Qr,xz \ I 1 A-^i^ A"" \ 
U,R - -Or^^^ 1 ^ij,23 = -A^- 1 A''^^' with g,, > 1, (3.41) 
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up to the first order of O^^ij = 1 ~ 3). After tlie seesaw meclianism, we obtain the 
Majorana mass matrix of light neutrinos in low-energy as, 
















1 


-1 





-1 


1 






Ai6 -Ai2 

A^ -A^ 1 

2A16 X\l - A^) 

-2A8 X\l + A-^) 

A8(l-A'^) + -2A4 

2 A^ + A^ A^-A^ 
A^ + A^ 2X^ 
A^ - A^ -2A^ 

2A8 A^-A^ A^ + l 
A^ - A^ -2A^ 1 - A^ 
A^ + 1 1 - A^ 2 



(3.42) 
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When we operate the Vtb to My as V^^MyVTQ, the neutrino mass matrix is 





24 











r')22 













1 










) 




. 






—\/3cic+ 






\/6C2C+ 


) 










C3A« 






2x^C2A^ 







-a/2CiC2 

— -\/3cic+ a/6c2C+ 

-2V2C1AS 

C3A« 

v^c_(A^ + A8) -v^c_(2A4- A^) 

V2X^ 
y2A8 2A8 -v^A^ 

-Vqx^ 

-ci ^ 
— ci 2a/2c2 \/3c_| 
v^c+ 



v^c_(A4 + A8) 
v^c_(2A'^- A«) 
-6v/2c+A^ 



(3.43) 



where 



ci = 1 - A^ - 2A« 



C2 



1 - A^ + A^ C3 = 1 - A^ + 4A« 



1±A^ 



(3.44) 



This mass matrix can be rewritten by 

M = Mo + Mos = Mo 



m 
m 



R 



12 



ml 



m 
m 



''13 ^23 '"'3 

where Aio comes from the diagonal elements of M^, which is given by 



R 



m 



R 
13 
R 
23 
R 



(3.45) 



\ 



-V2mi + ^^7713 

rhi + ^^^Thg, fh2 + 2mi + ^7713 
V^rhi - ^rh-s \/6mi - ^rha 



mi + ^ms 



(3.46) 



In ( 13.45P A^off has effects from the off-diagonal elements of M/j. In order to obtain ex- 
perimentally accepted mixing structure without unnatural cancellations, we focus only on 
a case that the collections from the off-diagonal elements of are small enough not 
to spoil the nearly tri-bimaximal mixing constructed by the cascade neutrino Dirac mass 
matrix. This means that the resultant structure of neutrino mass matrix given in fl3.43p 
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Xi 


X2 




Mr/M 


17 


10 




( A« A« A« \ 
A^ A^ -A^ 
A« -A^ 1 ) 






( X"' A" Ai" \ 
A^^ A^o A^ 
A^o A^ 1 ) 




18 


11 




( A« A« A« \ 
A8 A^ -A-i 
AS -A^ 1 ) 






( Ai« A^y All ( 
A19 All ^8 
All ^8 1 y 













Table 3: The textures of the neutrino Dirac mass matrix of the cascade form and the right- 
handed neutrino Majorana one of the non-diagonal from. The matrices are constrained by 
the experimentally observed values of the neutrino masses with the condition xi — X2 = 7. 

should not be drastically differed from the fl3.46p . and thus the magnitude of neutrino mass 
eigenvalues fl3.25p -( l3.27p and the above four constraints should be satisfied at the leading 
order even in non-diagonal Mr case. These discussions give the following neutrino mass 
eigenvalues up to the next leading order 



i 



mi ~ — ^ + mf = mi + mf , (3.47) 
6M 

ma ~ (^3Ai6""^ + ^ + mf = m2 + m^ + 2mi, (3.48) 

mg ~ (^2AS-^^ + 1^ ^ + mf = mg + mf + ?,rn^, (3.49) 
where mf include effects from the off-diagonal element of Mr described by 

mf ^ ^(A--^|23-2A«--^«,i3), (3.51) 

mf = ^A-^H2A%,i2 - ^i?,i3)'- (3.52) 

Typical textures of non-diagonal Mr are presented in Tabs. [3] and HI The collections to 
the generation mixing angles are also estimated as 

0<'' ^ (3^53) 



m2 + m2 



V^mi - ^m3 + -"^^ 



^ - W-'-'"^ (3,54, 

(m3 + mg J — (m2 + ma) 



mg + mg 



''Detailed discussions is given in the Appendix. 
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Xi 



17 



X2 II M^p/Vu 

A« A« A^ 
A^ A^ -A^ 
A« -A^ 1 



A^^ 


Ai« 


Ai" 


A^s 


A^ 


A^ 


A^o 


A6 


1 


Ai« 


Ai« 


A^" 



10 



A« A8 A^ 
A8 A^ -A"^ 
AS -A^ 1 



A18 
Aio 



Aio 
A^ 



A^ 
1 



Table 4: The textures of the neutrino Dirac mass matrix of the cascade form and the 
right-handed neutrino Majorana one of the diagonal from. The matrices constrained by 
the experimentally observed values of the neutrino masses with the condition xi — X2 = 8. 



where 



<2 - -^[^^'^^23^3- 2(2^,^,12 -A-%,i3)m2], (3.56) 



m 



R 



1 



-^[0R,23m3 + A^S(2A''^^«,i2 - ^i?,i3)(l - A-«eK,i3)m2], (3.57) 



A~^ 4 



(3.58) 



^ I ^(2 + X-^eR^23)0R,23rh3 - ^A-X,i2'^2 

Finally, the PMNS mixing angles including collections from the off-diagonal elements are 

mi[3(m3 + mf ) — (m2 -|- m^)] 



sin ^ 



23 



a/2 a/2 (m-s + mf ) [(ms + ) - (m2 + m^)] 

A'^^~'^2 m3(m2+"l^) 1 771^3 



3a/2 (?Ti3 + m|^) [(777,3 + TTT'i) - ("^2 + m2) a/6 ?Ti3 + mf 



sin 6 



13 



a/3 (m3 + ) - (777-2 + rn^) ' 

_)^di-d2 jTjg [(,7^3 + ,72^) _ I (^77^2 + 77i^)] 

y/2 {rhs + 7?7f ) [(7713 + 77if ) - (7712 + rn§)] 

■\/27?7i(7772 + TTlf ) 2 



(3.60) 



77i 



13 



(7713 + 777^) [(7713 + 7?7f ) - (7772 + 777^)] ^6 7773 + ^3 



777; 
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a/3 (77I3 + 777^) - (7772 + 77lf ) 

3.4 Charged lepton and quark sectors 



(3.61) 



At the end of this section, we study the charged lepton and quark sectors. Under the 
condition of the 5*0(10) scenario, we examine quantitative features of the masses and 
mixing angles. 
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We take the charged lepton mass matrix as 



Me ^ I 6d -3\d -3Xd I ^dVd, (3.62) 



-3A 



d 




in our study. The magnitudes of cascade parameters can be partially evaluated from 
the observed values of charged lepton masses, and are given by \Xd\ — m^/(3mT-) and 
\6d\ — Sy/m^ffT^/mr. We find that the corrections from the charged lepton sector are 
generally small; the total leptonic mixing angles can be written as 



sm6'i2 ~ ^ + ^— — ■ R^ + W , (3.63) 



. ^ 11 mi[3(m3 + m^) — (m2 + m^)] 

a/2 a/2 (ms + mf ) [(ms + ) - (m2 + m|^)] 
A'^^~'^2 rh^{rh2 + 711,2) 1 



3a/2 (ms + mf ) [(ma + ) - (m2 + m^) a/6 "^3 + "^f 

+ 7—^, (3.64) 



sin 6' 



13 



+ 



a/3 (ma + mf ) - (m2 + mf ) 3m^ 
\d1-d2 [(^^2 _^ _ l(^m2 + m^)] 

^^2 (-^3 + rnf) [("^3 + 'm§) - (m2 + mf )] 
-\/2mi(m2 + m^) 2 



mi3 



R 



(ma + mf ) [(ma + mf ) - (m2 + mf )] ^ ma + m; 
1 mS, 3 Irrip 



a/3 (ma + mf ) - (m2 + m^) a/2 y m 



at the first order of perturbations. 

Next, we comment on the quark sector. One must remember that the mass matrix of 
the H.C form is motivated for the mass spectra and mixing angles of quark sector. The 
mixing matrices for the up- and down-sector are given by the cascading mass matrices 
fl3.19p and f l3.20p . From the mass matrices, one can estimate the following mixing angles 



0{1) 0{\) 

Vd = I 0{\) 0{l) 0{\^) I , K = I OiX") 0{l) 0(A4) I , (3.66) 





C)(A3) 0{\^) 

where Vd and Vu are unitary mixing matrices determined by and Md- It can be easily 
seen from the structure of Vd that the experimentally observed values of CKM matrix 
can be realized at the leading order and the collections from the Vu are generally small. 
Detailed numerical calculations are given in the next section. 
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4 Phenomenologies 



In this section, we perform the numerical study of phenomenologies based on the above 
analyses of the cascade textures for the quark and lepton sectors: the PMNS mixing 
angles, lepton flavor violation (LFV), baryon asymmetry of the Universe (BAU) via thermal 
leptogenesis. 

4.1 PMNS mixing angles 

Firstly, we show numerical analyses of the generation mixing angles of the quark and lepton 
sectors predicted from the cascade model. Here, we investigate two typical types of minimal 
texture for the neutrino Dirac and right-handed Majorana neutrino mass matrices, 

(CyX^ CyX^ CyX^ \ 

c^A« KX^ -KX^ Vu, (4.1) 
cA^ -KX^ a, J 

/ fnX'' enX'' dnX'' \ 
Mr = I e^A^^ crAIO bRX' I M, (4.2) 

and 





cA' 


cA' 


cA' 


bA" 


-bA' 


cA' 


-bA' 


a„ 




crX'' 


rf^A^o 


CrX^' 


crX'' 


bRX' 


dRX'^ 


bRX' 




cA' 


cA' 


cA' 


cA' 


bA' 


-bA' 


cA' 


-bA" 






crX^' 


dRX''^ 


CrX'' 


crX'' 


bRX^ 


dRX'' 


bRX' 





Model II : M^n = cA^ bA^ -bA^ v^, (4.3) 

\ c^X^ -bA'^ au / 

f JrX'' crX'' dRX''\ 
Mr = crX^^ crX'^ M, (4.4) 

\dRX^^ aR J 

for the cases of the condition xi — X2 = 70 Here a,^, bi, c^, and qr, ■ ■ ■ , /r are complex 
numbers whose absolute values are taken as 0.4 ~ 1.4. In both models, the following 
charged lepton, up and down quark mass matrices are utilized 

CeA^ deX^ 
Me = I eeX^ -3ceA2 SbeX^ I XVd, (4.5) 



dA^ —SbpX^ ctp 



and 




Mu = cA"" bA^ -bA" Vu, Md = eA^ cA^ bA^ Xv^, (4.6) 




^We focus on only the xi — X2 — 7 case given in Tab. [3] as a typical example (Model I). In order to see 
effects from off-diagonal elements of right-handed neutrino mass matrix to the PMNS mixing angles, we 
also analyse a slightly different model for Mji as a comparison (Model II). 
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Figure 1: Predicted PMNS mixing angles from Model I and II: The solid and dashed lines, 
and white background correspond to the best fit value, boundaries of la errors and the 
region within 3a. 

where Og, ■ ■ ■ , Ce are also complex values whose range are the same as a,^, b^, Cu, and qr, 
• ■ ■ , fn. Moreover, notice that the mass matrices of down and up quarks have the same 
structures of ones of charged lepton, except for the GJ factor, and neutrino respectively 
because of the mass relations of 5*0(10) model at the GUT scale. 

The results of numerical calculation for the PMNS mixing angles in Model I and II 
are given in Figs. [1] The all plots can fit the data of quark masses, CKM mixing angles, 
charged lepton masses, and the mass ratio of two mass squared differences of neutrino. 
These are 

TTl Tfl 

— = 0.0026(6), — = 0.0023(2), yt = 0.51(2), (4.7) 
rric mt 

— = 0.051(7), — = 0.018(2), yfe = 0.34(3), (4.8) 
rris irib 

A = 0.73(3), Ac = 0.227(1), p = 0.22(6), 7^ = 0.33(4), (4.9) 

at the GUT scale, where the numbers in parentheses mean an uncertainty in last digit |16j|ii1 

I'Wc have taken a^, and a^, as 0.49 < < 0.53 and 0.31 < a^X < 0.37 to fit yt and yf, in our numerical 
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Here A, Ac, p, and fj are the Wolfenstein parameters which correspond to three mixing an- 
gles and one phase in the CKM matrix; yt and yb represent largest eigenvalues of Yukawa 
matrices for the up-type and down-type quarks, respectively. These values are calcu- 
lated with the 2-loop gauge coupling and 2-loop Yukawa coupling renormalization equa- 
tion taking tan/3, threshold corrections 'yt,b,d, and an effective SUSY scale as tan/3 = 38, 
7fe = —0.22, 7rf = —0.21, 7f = 0, and msusY = 500 GeV, respectively. The threshold 
corrections are approximated by 



otan f3 fiAt ^ 4 2tan/3/iM3 4 2tan/3/iM3 /, 

I a 



where /i, Af, mg, g^, and M3 indicate supersymmetric Higgs mass fiHuHd, soft top quark 
tri-linear coupling, mass of the squark q, strong coupling, and gaugino soft breaking mass, 
respectively [17]. The utilized values of SUSY can lead to the relations at the GUT scale, 

— = ^ = -^ = 1, (4.11) 

m-j- SrUe 

in a good accuracy. Therefore, the above values can automatically reproduce the experi- 
mental observed charged lepton masses at the low energy. 

The numerical calculations of Model I suggest that the predicted region of solar angle 
covers the experimental upper bound but the model gives a predicted lower bound around 
0.29 < sin^^i2. On the other hand, a constrained region is predicted for the reactor angle 
as 0.002 < sin^^i3 < 0.007. This result is one of important predictions of the present 
cascade textures. It might be checked by the upcoming DoubleChooz [IB], RENO [TU] . 
and DayaBay [20] experiments as the reactor experiments in addition to the accelerator 
experiments such as T2K [21] and NOi^A [22] o Finally, the atmospheric angle covers 
the current experimentally allowed region. However, there is a relatively clear correlation 
between the magnitudes of reactor and atmospheric angles in Fig. [H On the other hand, 
the result for Model II can cover the experimental allowed region because of largeness of 
corrections from the right-handed neutrino mass matrix but there is an upper bound of 
sin^ 6'i3 which is sin^ ^13 < 0.015. Therefore, we can conclude that the minimal cascade 
textures in the context of SUSY 5*0(10) (Model I) can lead to clear predictions, which 
are 0.29 < sin^ 9i2 and 0.002 < sin^ 6^13 < 0.007, and relatively sharp correlations between 
the reactor and atmospheric angles. In the next-to minimal cascade model (Model II) lead 
to only the upper bound of the reactor angle while the model can explain about almost 
ranges of PMNS mixing angles. The minimal realization of cascade model is predictive 
and interesting in the framework with the cascade hierarchies in SUSY 50(10) GUT, and 



analyses since the values of yt and yb are almost determined by ones of a,/ and Og, respectively. 
**See j23| for an excellent review of sensitivities of the upcoming experiments. 
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thus it might be checked by the future experiments. It would be too difficult to distinguish 
other cascade textures even if the future precision data of neutrino oscillation experiment 
could be used. It is worth studying a new method to check the models. 

At the end of this subsection, we give a brief comparison between our results and ones 
from a similar 5*0(10) approach, which utilizes type II seesaw mechanism and a simple 
ansatz such that the dominant Yukawa matrix has rank one [2l]|ll| The work gives some 
predicted regions for leptonic mixing angles based on three typical models in this approach: 
(A) Vy = 1 case where Vpmns = VeVrBVj, and Vg and are diagonalizing matrices of 
Yukawa matrix for the charged lepton and neutrino mass matrix for light neutrinos in a 
tetrahedral coordinate, respectively, (Bl) K ^ 1 and (/*^*''')i2 = {/'"''^i = (f ''*'^)i3 = 
^ytetra-j^^ _ g ^]-^grg jtctra jg ^ coupliug to 126+126 Higgs in tetrahedral coodinate, and (B2) 

K ^ 1 and (/*^*''^)i2 = (/*"*''^)2i = (/*'*''')23 = (/*'*'^)32 = 0. The numerical calculations 
predict sin^ Ou ~ 0.28 and 0.006 < sin^ < 0.008 in the model A, 0.32 < sin^ 612 < 0.33, 
0.003 < sin^^ig < 0.006, and 0.59 < sin^ ^23 < 0.60 in the model Bl, and sin^ ^12 ~ 0.30, 
0.04 < sin^ < 0.06, and 0.62 < sin^ 623 < 0.63 in the model B2 by using the data 
Am2i/Am|;^ = 0.027 — 0.038 at 2a level. These can be compared with our results and 
might be also checked by the future neutrino experiments. 

4.2 Lepton flavor violation 

Next, we investigate the branching ratios of LFV process k — )■ Ij'-f, in our cascade lepton 
mass matrices. We suppose that soft SUSY breaking masses of sleptons are universal at 
the GUT scale, Agut? for simplicity. In the case, the off-diagonal matrix elements are 
generated by radiative corrections from the Yukawa couplings of neutrinos [25]. The one- 
loop renormalization group evolution gives the left-handed slepton masses. The leading 
contribution is estimated by 



-.^< + \^<^\\j2{Ml^MM.o),A^(]^) (for.T^j), (4.12) 



H Jij 



Sir'^v'^ sin^ f3 



where mo and ao are the universal SUSY breaking mass and three-point coupling of scalar 
superpartners given at the GUT scale. The branching fractions of each LFV process are 
roughly given by 



B,.(,.^,,,):.|^Mk«:tan=/?, (4.13) 



^SUSY 



t^The paper by Dutta, et al. in refs. |3] presents an 54 flavor model, which is one of rcahzations of the 
rank I approach by flavor symmetry. Here we focus on the general results of this approach given in [24] . 
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Br(/i e-i)/B 


Br(r e-i)/B 


Br(r ^ /i7)/5 


Ml [GeV] 


Ms [GeV] 


Mg [GeV] 


6.95 X 10^19 


1.57 X 10^1^ 


2.23 X 10^12 


3 X 10^ 


8 X 10^ 


2 X 10^6 



Table 5: Typical magnitudes of branching ratios for lepton flavor violating rare decay 
process. 

where a, Mw, and msusY are the fine structure constant, W boson mass, and a typical 
mass scale of superparticles, respectively. These branching ratios are estimated by 

2 



Br(/x — )■ 67) 



3a 



B 



A^^ln 



Br(r — i- 67) 



Br(r — i- /i7) 



3a 



B 



3a 



-5 



In 



|Mi| 
Agut 

|Mi| 
Agut 

iMil 



A^^^ln 



A 



GUT 



A«ln 



IM2I 
Agut 

IM2I 
Agut 

IM2 



A^^ln 



A 



gut 



A«ln 



A^n 



IM3I 
Agut 

IM3I 
Agut 

IM3 



n 2 



(4.14) 

(4.15) 
(4.16) 

kol = 

'^susY- Typical magnitudes of the branching ratios are shown in Tab. El In these analyses, 
Agut = 2 x 10^^ GeV is taken. These results are compared with the current experimental 
upper bounds at 90% confidence level [26], 



A, 



gut 



for both Model I and II. Here we define B = (Mty/?7isusY)^fan^/3 and take rriQ 



Br(/i ^ 67) < 1.2 X 10"", Br(r 67) < 1.2 X 10~^ Br(r ^ /i7) < 4.5 x 10^^(4.17) 

The magnitudes of the branching ratios for the lepton flavor violating process in the model 
with the applicable heaviest right-handed Majorana neutrino mass are far below the exper- 
imental limit. Once one fix the value of tan/3, the current experimental limit gives lower 
bound on msusY, which is e.g. msusY > 41.6 GeV from Br(r — fi'-f) < 4.5 x 10~^ with 
tan/3 = 38. In the present case, the bound can be easily satisfied. 

4.3 Leptogenesis 

At the end of this section, we examine whether the thermal leptogenesis [28] works in our 
model. The CP asymmetry parameter in the decay process of right-handed neutrino, Ri, 
is given by 



E, r(i?. ^ L,H) + E, nR^ ^ ^Hr) ' 



(4.18) 



where Li and H denote the left-handed lepton and Higgs fields. An approximation for 
at low temperature is estimated by [22], ei = ^ Ei^i I'^[(^ii)^]-^('^i)/l^ii|) where = 
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|Mi/Mip, A = {DM^oMI^D"^) /v^, and D being a diagonal phase matrix to make the 
eigenvalues Mj real and positive. The function F describes contributions from the one- 
loop vertex and self-energy corrections, 



We here define the resultant CP asymmetry, ?7cp, as the ratio of the lepton asymmetry to 
the photon number density n^, rjcp = 135C(3)Ksei/(47r^(7^,n^). In the equation, n, s, and 
are the efficiency factor, entropy density, and the effective number of degrees of freedom 
in thermal equilibrium. They are given by [30], s = 7.04n^, = 228.75, and 



BAU, rjB, is transferred via spharelon interactions as tjb = ~8?7cp/23. Finally, the baryon 
asymmetry in our model is predicted as ?7b ~ 10^^^ sin where 6b = 63 — 6i and 6i = 
arg(Mj). These results are compared with the current observational data at 68% confidence 
level from the WMAP 7-year resulting in the standard ACDM model [31j. We can see that 
the baryon asymmetry generated through the leptogenesis is too small to explain the BAU. 
This is because the hierarchy in the Dirac neutrino mass matrix in the 5*0(10) model is 
determined by the up-type quark mass matrix. Therefore, there is no freedom to adjust 
the Dirac neutrino mass matrix such that the BAU can be generated in the present model. 
It is expected that enough BAU might be realized by extending our cascade model to the 
inverse seesaw case, see [32]. In the case, the structure of effective light neutrino mass 
matrix is slightly changed but the realistic PMNS mixing angels would be obtained. 

5 Discussion 

At the end of the paper, we give a comment on phenomenological aspects of proton decay. 
In general for SUSY GUTs, there are three sources that mediate the proton decays. The 
first one comes from the dimension-6 operators, arising from the exchanging of the heavy 
gauge bosons. Note that this type of operators exists in both non-SUSY and SUSY GUTs. 
These operators are significantly suppressed by 1/ Aq^t? therefore, we have no problem with 
the proton decay from these operators if Agut is large enough, i.e. Aqut > 10^^ GeV. The 
second source is from the dimension-5 operators, arising from the exchanging of color triplet 
Higgsino fields. In this case, the proton decay contribution is suppressed by 1/Mh, where 
Mb is the mass of the color triplet Higgsinos. To suppress the proton decay contributions, 
the mass Mh has to be very heavy, which can achieved by some doublet-triplet splitting 




(4.19) 
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mechanism [121 [131 113] • The third contribution arises from the dimension-4 operators, 
which are not suppressed by the GUT scale, however, these operators are ehminated by 
the -R— parity. In the class of 50(10) models, which do not contain the spinor 16 or 16 as 
the Higgs fields, the i?— parity is conserved. Otherwise, in order to avoid the proton decay 
contributions from these operators, the i?— parity has to be introduced by hands. Since 
we consider an SO{10) model without using the spinor Higgs representations, we have no 
problem with these dimension-4 operators [331 133] • 

6 Summary 

We have done texture analyses of cascade model in supersymmetric 5*0(10) model. The 
neutrino Dirac mass matrix of a cascade form can realize the tri-bimaximal mixing at the 
leading order while the down-quark one of a H.C form can lead to realistic structure of 
CKM mixing. This fact gives us a strong motivation to study cascade hierarchical textures 
in a grand unified theory. 

We analytically clarified possible structures of neutrino Dirac, charged lepton, quark, 
and right-handed neutrino mass matrices by estimating collection from them to the tri- 
bimaximal mixing. The numerical analyses based on two typical models have been also 
presented. The minimal cascade texture in the context of SUSY 50(10) GUT can lead 
to clear predictions for the PMNS mixing angles, which are 0.29 < sin^ 9i2 and 0.002 < 
sin^ ^^13 < 0.007, and relatively sharp correlations between the reactor and atmospheric 
angles. This result is a hot topic for the upcoming experiments of the reactor neutrino 
mixing angle. It might be checked by such future experiments. We have also shown 
that our typical cascade models can pass the constraints from the lepton flavor violation 
searches. For generating the BAU, we cannot generate enough asymmetry through the 
thermal leptogenesis mechanism in our model. Therefore, we need other mechanisms to 
generate the BAU. 
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A Constraints on structure of non-diagonal Mji case 

The constraints on the structure of non-diagonal M/j are presented in this Appendix. We 
defined the diagonalized mass matrix of the right-handed neutrino in fl3.40p and an unitary 
matrix which diagonalizes the Mr. The neutrino mass after the seesaw mechanism and 
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operating the Vtb is given in f l3.42p . All matrix elements are given by 
— \ 



2 1 \16-X2 n\8 

M22 ^ [3A''-^^ + - + + X-^'i-—eR,23 + el,,) + A-^M-2A«^^,i3 + ^^,13)], 

(A.2) 

A^33 ^ ^[2A«-^^ + ^ + A-^^(2AX,23 + ^)-2A^-^^efi,i2eH,i3 

+ A--(2A«^^,,2 + ^^U)]' (A-3) 
-M12 ^ --^[l + A«-^^^^,23 + 3A-^H-2A'%,i2 + A%,i3)], (A.4) 



3y2M 



1;! 



" ^^[l-2Ai2-^^ + 2A4-^^^^,23 



-^13 ^ -^-|^[l + A-^^(4Ai2 + 2A%,23 + 

- A--(2A^^H,i2^«,i3 - 4A^2^^,i2 - Ol,,)]. (A.6) 



+ 3A-^H2Ai2^^,,i2 - X'On^rs - 2A^^«,i3efl,i2 + ^l.is)], (A.5) 

2 



We require that the magnitudes of leading order of each term in this mass matrix are the 
same one as in the case of diagonal M/j case because the tri-bimaximal mixing can be 
already realized at the leading order by neutrino Dirac mass matrix. It leads to constraints 
on the mixing angles as follows: 

9R,n < -A^ -A-«+^\ 2A^+(^^-^^)/2^ ^A"^/^ (A.7) 

2 3 y/^ 

9r,23 < A^ ^A-^+^^ A^^/2 (A.8) 

Oru < ^A^12+x,^ ^(xi-X2)/2^ ^-6+x,/2^ 

6 2 

^/?.i2^i?,i3 < A^+-•-^^ ^A-^+^^ (A.IO) 

o 

After fixing the values of (rfi, ^2, a^i, 2:2) so that they must satisfy the four conditions, one 
can obtain the structure leading to maximal collections to the PMNS mixing angles and 
neutrino mass spectra as shown in Tabs. [3] and HI 
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